The quantization of the B = 1 and B = 2 Skyrmions 
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I. INTRODUCTION 
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We propose to set the Skyrme parameters F-n and e such that they reproduce the physical masses 
of the nucleon and the deuteron. We allow deformation using an axially symmetric solution and 
simulated annealing to minimize the total energy for the B = 1 nucleon and B = 2 deuteron. First 
we find that axial deformations are responsible for a significant reduction (factor of « 4) of the 
rotational energy but also that it is not possible to get a common set of parameters and e which 
would fit both nucleon and deuteron masses simultaneously at least for = 138 MeV, 345 MeV 
. and 500 MeV. This suggests that either m,r > 500 MeV or additional terms must be added to the 

■ Skyrme Lagrangian. 
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The Skyrme model []| is a nonlinear theory of pions that admits topological solitons solutions called Skyrmions. 
' '' These solutions fall into sectors characterized by an integer- valued topological invariant B. In its quantized version, 
O ^1 a Skyrmion of topological charge B may be identified with a nucleus with baryon number B. 

Since the model is non-renormalizable, a canonical quantization is not possible and one has to resort to semiclassical 
; quantization of the zero modes of the Skyrmion. This method adds kinematical terms (rotational, vibrational^.) to 
(— I ' the total energy of the Skyrmion. The B — 1 Skyrmion was first quantized by Adkins, Nappi and Witten 0, Q. It 
provided then a useful mean to set the parameters of the Skyrme model and e by fitting to the proton and delta 
masses. The experimental value of the pion mass ttItt completes the set of input parameters when a pion mass term is 
added to the Skyrme Lagrangian. The same calibration was assumed by Braaten and Carson 4] in their quantization 
of the B = 2 Skyrmions and their predictions of a rather tightly bound and small sized deuteron. Further analysis 
ly-j by Leese, Manton and Schroers 5] considering the separation of two single Skyrmions in the most attractive channel 
I led to more accurate predictions for the deuteron. 

' Yet all these calculations have three caveats: First, they used a rigid-body quantization, i.e. the kinematic term 
, is calculated from the solution which minimizes the static energy neglecting any deformation that could originate 

■ form the kinematical terms. This was pointed out by several authors p, 0| who proposed to improve the solutions by 
00 ! allowing the B = I Skyrmion to deform within the spherically symmetric ansatz. Second, even with such deformations 

■ it was noted that the set of parameters mentioned above does not allow for a stable spinning solution for the delta 
^ I since the rotational frequency fl would not satisfy the constraint for stability against pion emission, fi^ < |mj. These 

• ^ two problems were addressed recently in [1, Q . Assuming an axial symmetry, the calculations were performed using 
' a simulated annealing algorithm allowing for the minimization of the total energy (static and kinetic). It was also 
shown that the stable spinning nucleon and delta masses could be obtained only if the pion mass is fixed at more than 
twice its experimental value. One may argue that in this case m-^ could be interpreted as a renormalized pion mass 
which could explain its departure form the experimental value. A third difficulty remains: Fixing the parameters of 
the Skyrme model still involve the delta which is interpreted as a stable spinning Skyrmion although physically it 
is an unstable resonance. Recently, Manton and Wood [l3| took a different approach and chose data from B = 6 
Lithium-6 to set the Skyrme parameters. The purpose was to provide for a better description of higher B solutions. 
Unfortunately, their calculations were based on a two-fold approximation, the rational map ansatz and rigid-body 
quantization for values of B up to 8. 

In this work, we propose to set the Skyrme parameters such that they reproduce the physical masses of the nucleon 
and the deuteron. We allow deformation within an axially symmetric solution approximation and use simulated 
annealing to minimize the total energy for the B = 1 nucleon and B = 2 deuteron. We argue that this procedure 
provide a solution which is very close to if not the exact solution. Following the procedure in we find the sets 
of parameters Fj^ and e that are required to fit the nucleon and deuteron masses respectively for TOtt = 138 MeV, 
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345 MeV and 500 MeV. The numerical calculations are compared to those obtained with the rational maps ansatz 
and rigid-body quantization. We find the latter approximation to be misleading since it suggests that it is possible 
to simultaneously fit the nucleon and deuteron masses which is not the case when we perform an our numerical 
calculation even for larger ttItt = 500 MeV. However for the solution to remains realistic with regard with the size of 
the nucleon or deuteron, lower values of F^^ and e are to be excluded. 

In section II, we introduce briefly the Skyrme model and find the static energy for the axially symmetric solution 
proposed in The quantization of rotational and isorotational excitation using this solution leads to an expression 
for the kinetic energy terms at the end of section III. These expression suggest that the axial symmetry could be 
preserved to a large extent. Finally we discuss and compare our numerical results from simulated annealing with an 
axial solution on a two dimensional grid with that coming from rational maps with rigid-body approximation in the 
last section. 



II. THE SKYRME MODEL 



The SU{2) Skyrme Lagrangian density is 

p2 1 2 p2 

16 ^ ^ ' 32e 



Cs - Tr (L.Ln + -— Tr [L„ L.f + (TrC/ - 2) (1) 



where is the left-handed chiral current = U'^d^U and the parameters and e are respectively the pion decay 
constant and the dimensionless Skyrme constant. C/ is a SU{2) field associated to the pion field tt by 

U = (7 + iT - TV (2) 

where t are the Pauli matrices, tt = (tti, 7r2, tts) is the triplet of pion fields and the scalar meson field a satisfy 
cr^ -|- TT • TT = 1. The third term where m^r is the pion mass was first added by Adkins and Nappi [3] to account for 
the chiral symmetry breaking observed in nature. 

The field configurations that satisfy the boundary condition 

f/(r,t)^l for |r|^oo (3) 

fall into topological sectors labelled by a topological invariant 

^ = i/ '^'^d^ti^^} = '§^J d'^Tr (L,[L,,L,.]) (4) 

taking integral values. Skyrme interpreted this topological invariant as the baryon number. The minimal static 
energy Skyrmion for B — 1 and B — 2 turns out to have spherical and axial symmetry respectively. Since we arc only 
interested by these values of B, the solution will be cast in terms of cylindrical coordinates (p, 0, z) in the form 

a = ijj^ TTi = ^1 cos nO 112 = ^1 sin nO tts = 7/12 (5) 

introduced in |llj| where ip{p, z) = {ipi, ■02, "03) is a three-component unit vector that is independent of the azimuthal 
angle 9. The boundary conditions Q implies that — > (0, 0, 1) as + z"^ ^ 00. Moreover, we must impose that 
tpi = and dptlj2 — dpi^z = at p = 0. The B = 1 hedgehog solution appears as a special case of (O having spherical 
symmetry and corresponds to 

(01, ■02, V's) = (sin F sin 0, sin COS 0, cos F) 

where F — F{r) is the profile or chiral angle. 

With the axial ansatz ([5]) and a appropriate scaling ^ the expressions for the static energy and the baryon number 
become 

En = - f d^xCs 

(6) 



^ We have used 2-v/2/eF^ and Ft^/2\/2 as units of length and energy respectively. 
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B = - I dzdpiJi\dpi> X d,tp\ (7) 



with /3 = 2V|nzL. 

However, to describe baryons, Skyrmions must acquire a well defined spin and isospin state. This is possible only 
upon proper quantization of the Skyrmions as it will be done in the next section. As we shall see in the next section 
adding (iso-)rotational energy will in general brake the axial symmetry manifest in _B = 1 and B = 2 static solutions. 



III. QUANTIZATION 

Since the Skyrme Lagrangian |T]) is invariant under rotation and isorotation ^, the usual method of Skyrmion 
quantization consist of allowing the zero modes to depend on time and then quantize the resulting dynamical system 
according to standard semiclassical methods. From this perspective, the dynamical ansatz is assumed to be 

U{r,t) ^ A,{t)U{R{A2{t)r)A\{t) (8) 

where Ai, Ai are SU{2) matrices and i?y (A2) = 5Tr(ri A2Tj A\) is the associated SO{3) rotation matrix. Introducing 
this ansatz into the Skyrme Lagrangian ([T]) one gets the kinematical contribution to the total energy which can be 
cast in the form 

T = i a; Uij Qj - ai Wij bj + ^bi Vij bj (9) 
where aj — — iTr ^Tj AJAi^ , bj = ?Tr [tj A2A^^ and the C/y, Vij and Wij are inertia tensors 

Urj = - J d'a^jTr (T.T,) + ^Tr ([i,., T,] [L,, (10) 

Vij = - (^^y-^ <^ikl£]mn J d^XXkXm^Tl (LiLn) + ^Tr([Lp,L;] [Lp,L„])|, (11) 



W,, = (^\e,ki J dlTa;fc|Tr(r,iO + ^Tr {[L,n,T,] [L„,,L„])| 



(12) 



where Ti = ill'' [^-jU]. For the axial ansatz ([5]), these tensors are all diagonal and satisfy Uu = U22, Vu ~ V22, 
Wii — W22 and U33 — — The components of these inertia tensors are 

Uii = 2tt (^^^ J dzdppi^i,l + 2^1 + i {dpiP ■ dptP + d.iP ■ d.xP + n^^^ + {dp^sf + [dz^sf 

ii>V' ^^^^ 
P 




U33 = 2n / dzdppyjf {dp^ ■ dpxP + 8,1^ ■ d.i^ + 2) , (14) 



^ Since we are interested only in the computation of the static properties, we will ignore translational modes and quantize the Skyrmions 
in their rest frame. 
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- zd,^\-^ [ 1 + n'^^ + z'n'^ + ^ (p' + I^p^ x ^^^I'}. (15) 



'4^3 



{dzi^zf + {dp^^f + % 



(16) 



2p 



(2 + dpxP ■ dpxp + dzip ■ dzip) 



Let us note that Wn ^ only for n = 1. In order to obtain the energy corresponding to the nucleons and the 
deuteron, we must compute the Hamiltonian for the rotational and isorotational degrees of freedom in term of the 
inertia tensors (fT3|) to (fT6|) as well as its eigenvalues for the states corresponding to these particles. 
The body-fixed isospin and angular momentum canonically conjugate to a and b are respectively 



dT 

Ki = — ^ U^Ja■J - W^JbJ, 



(17) 



_ dT 



-W^,a,+V^jbj. 



(18) 



These momenta are related to the usual space-fixed isospin (I) and spin (J) by the orthogonal transformations 

= -R{Ai),,Kj, (19) 



J, = -R{A2)lL,. (20) 
We can now write the Hamiltonian for the rotational and isorotational degrees of freedom as 

= K • a L • b - T 

where T is the kinematical energy ([9]). 

The quantization procedure consist of promoting four-momenta as quantum operators that satisfy each one the 
SU{2) commutation relations. According to and ([^0]) . we see that the Casimir invariants satisfy — and 

—li^. Then the operators form a 0(4)i_k ^ 0(4)l,j Lie algebra. The physical states on which these operators act 
are the states formed in the base \^) = lH^h^lljjsh) where —i < 13, —j < ja and ^3 < j that satisfy the constraints 
formulated by Finkelstein and Rubinstein [12| . One of these constraints, namely, 



g2^»„.L|^^ = e^'^'n-Kl^^ = (-1)^1^'), (21) 

implies that the spin and isospin must be an integer for even B or a half-integer for odd baryon numbers. Another 
constraint 

(nXg + ^3) |*> = (22) 

comes from the axial symmetry imposed on the solution here. 

This formalism allows one to obtain the total energy in terms of the inertia tensors of (|13m6p for the nucleon 
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and for the deuteron [J| 

Ed = + (24) 

Vll 

For the B = 1 solution the minimization of the static energy Ei leads to a spherically symmetric solution. The 
nucleon mass in (j23p however receives (iso-) rotational energy contributions from (iso-)rotations around the three 
principal axis among which two are equivalent (direction 1 and 2). The deformation should then lead to a axially 
symmetric solution along the z-axis (direction 3) as argued in Refs.[l, 0. 

The situation is somewhat different for B — 2 which is known to have a toroidal solution upon minimization of the 
static energy i?2- In our scheme it corresponds to the axially symmetric solution with respect to the z-axis. Assuming 
such a solution the deuteron mass in (|24[) gets kinetic energy contributions which may be recast in the form 

^™' = 2^+2^=^l^' (^^^ 

which means that the contributions come from rotations perpendicular to the axis of symmetry. Unfortunately there 
is no guaranty that in a non-rigid rotator approach the axial symmetry will be preserved when the rotational terms 
are added. However the deformations are not expected to be very large since the magnitude rotational energy only 
accounts for less than 4% of the total mass of the deuteron in the rigid-body approximation |4(] . Allowing deformations 
should increase the moments of inertia and bring the relative contribution of rotational energy to an even lower value. 
Clearly large changes in the configuration are prohibited by an eventual increase in the static energy E2 so, from that 
argument alone we can infer that deformations from axial symmetry are bound to be contained into a 4% effect. Our 
results will show that non- axial deformations represents less than 1% of the deuteron mass if they contribute at all. 
This justifies the use of the axial solution in ([5]). 

In order to obtain the configurations ij}{p,z) ~ (-01, ■02, "03) that minimize the energies defined in (|23p and (|24p . 
we use a two-dimensional version of a simulated annealing algorithm used in f\3\. The minimization is carried out 
on a grid in a plane (p, z) made up of 250 x 500 points with a spacing of 0.042 ^. The algorithm starts with an 
initial configuration ■j/'o(yO, z) on the grid and evolves towards the exact solution. Here ■j/'o(p, z) is generated using the 
suitable rational map ansatz to ensure that the initial solution have the appropriate baryon number, i? = 1 or 2, 
with a profile function of the form 

F(r) = 4arctan(e""'') (26) 

as inspired from Ref . [l5l| . Here a is a parameter chosen so that the entire baryon number density (such that B = 1 
or 2 respectively for the nucleon or the deuteron) fits into the grid. 



IV. RESULTS AND DISCUSSION 



There are several ways to fix the parameters of the Skyrme model e, F-^ and m-^ and the authors of Ref. [Sj] showed 
that one must take some precautions in order that the spinning solutions for the nucleon and the delta remain stable 
against pion emission. And indeed it was found that, in order to achieve a fit for the energies of spinning Skyrmions 
to the masses of the nucleon and delta, one must impose a value for the pion mass that is larger than its experimental 
value. Differences between the fitted and experimental values of and ttItt are not proscribed since after all, the 
values that enters the Lagrangian ([1]) are the unrenormalized parameters which could differ from the physical ones. 

Even so it remains that this procedure still assumes that a physically unstable particle, the delta resonance, is 
described as a stable spinning Skyrmion and to be perfectly consistent one should instead rely on stable particles. For 
these reasons we chose to carry out calculations using data from two stable particles, the nucleon and the deuteron. 
We proceed as follows: Assuming a value for the pion mass and e, and an initial value for the Skyrme parameter 
Ft, we compute the lowest energy solution for the spinning Skyrmion corresponding to the nucleon using simulated 
annealing which leads to a mass prediction. We iterate the procedure adjusting values of -FV until the predicted mass 
fits that of the nucleon. The same procedure is repeated for the B = 2 deuteron as well as for several values of e. A 
set of points requires about two weeks of computer calculations on a regular PC. Performing an equivalent calculation 
on a 3D grid with similar spacing for a non-axial solution for example would be prohibitive. This explain in parts 
why we use the axially symmetric ansatz. 



^ These parameters were adopted in order to be similar with those used in 
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Initially, the first set of calculations was performed with the pion mass equal to its experimental value m,r = 138 
MeV. Since Ref. [1| suggests that a pion mass value of = 345 MeV or more is necessary to avoid instability due to 
pion emission we repeated the calculations and evaluated the Skyrme parameters using TOtt — 345 MeV. The results 
arc displayed in FIG. [TJ Although it may seem interesting here to consider values of e larger than those illustrated 
on FIG. [U the physical relevance of lower values of e is questionable. Below a certain value, the rotational energy of 
the Skyrmion is larger than the contribution coming from the pion mass term. As was highlighted in Refs. [1, [3|j this 
leads to an unstable Skyrmion with respect to emission of pions. 




FIG. 1: Ftt as a function of e for which M + E^ot is equal to the nucleon mass (circles) and that of the deuteron (squares). 
The solid and dashed lines correspond to the results obtained from the rigid body approach for the nucleon and deuteron 
respectively. The set of data and lines at the upper left are for m,r = 138 MeV whereas that at the bottom right are for 
rUn = 345 MeV. 

Our results for B = 1 reproduce the same behavior as in Refs. [1, Q, i.e. the deformation of the Skyrmion becomes 
relatively important only for larger values of TOtt and e. However, the results for B = 2 are far more interesting 
with respect to deformation. Indeed, as we can see directly from FIG.[T1 there is a noticeable difference between our 
numerical results (squares) which allow for deformations as long as they preserve axial symmetry and that of the 
rigid body approximation (dashed lines). Note that our implementation of the rigid body approximation shown here 
relies on the rational map ansatz which is known to be accurate to a few percent and so it could not be responsible 
alone such large difference in energy. This difference indeed corresponds to a surprisingly much smaller energy for the 
deformed spinning Skyrmion than what is obtained from a the rigid body approximation wether it is based on the 
rational map ansatz or not. 

To illustrate this difference, we carried out our minimization procedure using the set of parameters -F^ — 108 MeV 
and e = 4.84 for the deformed deuteron. The results, set SA, are listed in of Table HI For comparison we also 
present three rigid body calculations: set SA-RB performed with our simulated annealing algorithm, set Ref.Q-RB 
from Braaten et al and finally set RM-RB obtained through the rational maps approximation. Note that our results 
SA-RB agrees fairly well with that of Ref.[3|. Clearly the mass of the deuteron En must be larger than iJ™'" = 1655 





SA 


SA-RB 


Ref.f4]-RB 


RM-RB 


Exp. 


Ed (MeV) 


1679 


1716 


1720 


1750 


1876 


E,ot (MeV) 


13.3 


60.6 


61.2 


55.2 




{r')T (fm) 


0.94 


0.93 


0.92 


0.94 


2.095 



TABLE I: Total energy, rotational energy and charge radius of the deuteron using parameters F,r ~ 108 MeV and e = 4.84. 
The results are presented for our simulated annealing calculations with axial deformation (SA) along with three rigid body 
calculations: with simulated annealing (SA-RB), from Braaten et al. (Ref.Q-RB) and with the rational maps approximation 
(RM-RB) respectively. The last column (Exp.) shows the experimental values where the charge radius come from [161. 

MeV, the minimum static energy for this choice of parameters. Computing the rotational energy for this solution leads 
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to the rigid body approximation result (SA-RB) for the deuteron mass of Eu = 1716 MeV. As expected the results 
for the deformed deuteron (SA) Ejj = 1679 MeV lies between these two values. Moreover our axial solution brings 
the relative contribution of the rotational term to about 1% of the minimum static energy ii^™™. So allowing for axial 
deformation reduces by about a factor of 4 the rotational term going from 60.6 MeV to 13.3 MeV. As for non-axial 
deformations that might be present in a completely general solution, they must at the very most represent an 1% 
correction to the deuteron mass since it is bounded by E^^™ (= 1655 MeV) < £;g^act < -^gjg ]y[ev). On the 

other hand they may still represent a significant portion of the remaining 13.3 MeV rotational energy. We conclude 
nonetheless that the axial symmetry ansatz represents a very good approximation of the deuteron configuration. 
In addition, since these bounds are both based on axially symmetric solutions and their energy only differ by 1%, 
one can even contemplate the possibility that the exact solution may have axial symmetry and therefore is to our 
solution, contrary to what physical intuition might suggest. This has yet to be proven and needless to say that such 
a demonstration would require a 3D calculation with a level of precision less than 1%. 

Note that, both for = 138 MeV and 345 MeV, the solid and dashed lines intersect for a relatively low value of 
e. However despite what this rigid body calculations suggests, there is no set F-^ and e that leads simultaneously to 
the masses of the nucleon and the deuteron. However, the data of FIG. [1] would indicate that the gap between the 
values of -FV decreases for a larger pion mass. For that reason, we repeated the calculations with a relatively large 
TOtt = 500 MeV. Some of the computations required that we adjust the spacing of our grid to 0.0057 for smaller e 
since the configurations turned out to be much smaller in size. The results for to^ — 500 MeV are presented in FIG. 
[5) The main conclusion one can draw is that even for this pion mass, not common set of Skyrme parameters can be 
found. 

It is also interesting to note that the gap between the values of decrease as e decreases. This would suggest that 
at very low values of e a fit is possible. However computing the charge radius, i.e. the square root of 



leads to an significant increase for the radius for small values of e and F^r as illustrated in FIGS. [3]to[5l So lower 
values of e are incompatible with the physical size of the deuteron and nucleon and may be discarded. The results 
also indicate that the best fit for the radius of the nucleon and deuteron would favor intermediate values of e around 
e ~ 3.5 while it looks fairly insensitive to large changes in m^. 

To summarize, our calculations showed that the axial symmetry ansatz is a very good approximation of the exact 
solution for the deuteron. This hints at the possibility that it may even represent the exact solution. This remains 
to be proved with a general 3D calculation. We also found that allowing for axial deformation reduces the rotational 
energy by a significant factor. On the other hand we found that it is not possible to get a common set of parameters 
and e which would fit both nucleon and deuteron masses simultaneously at least for TOtt — 138 MeV, 345 MeV and 
500 MeV. This conclusion should hold even for the exact B = 2 solution since if the solution was allowed to adjust 
free of any symmetry contraints it would achieve a configuration with lower total energy which would require larger 
values of for the same set of e. This suggests that either ttIt,- > 500 MeV or additional terms must be added to the 
Skyrme Lagrangian. We also observed an increase in the deformations due to the spinning of the B = 2 Skyrmion 
(deuteron) especially for larger values of e and -FV so the rigid body approximation may not be appropriate in that 
case. 

This work was were supported by the National Science and Engineering Research Council. 




(27) 
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FIG. 3: Charge radius of the nucleon (bold circles) and deuteron (bold squares) for the set of parameters of FIG.[T]with m,r = 138 
MeV. The curves and the horizontal lines, solid for the nucleons and dashed for the deuteron, correspond respectively to the 
results obtained from the rigid body approach and to the experimental data. 
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FIG. 4: Same as FIG. E] for = 345 MeV. 




FIG. 5: Same as FIG. E] for = 500 MeV. 
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